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We study the entanglement generation of operators whose statistical properties approach those 
of random matrices but are restricted in some way. These include interpolating ensemble matrices, 
where the interval of the independent random parameters are restricted, pseudo-random operators, 
where there are far fewer random parameters than required for random matrices, and quantum 
chaotic evolution. Restricting randomness in different ways allows us to probe connections between 
entanglement and randomness. We comment on which properties affect entanglement generation 
and discuss ways of efficiently producing random states on a quantum computer. 
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Highly entangled, random, quantum states play a cen- 
tral role in many aspects of quantum information process- 
ing (QIP). Protocols enabled by random quantum states 
include superdense coding remote state preparation 
[2j , data hiding schemes [3j , and single spin measurement 
4] . Random states are produced from computational ba- 
sis states by applying random unitary operators. How- 
ever, the implementation of operators randomly drawn 
from the circular unitary ensemble (CUE), the space of 
all unitaries, is inefficient. 

Independent of its usefulness in QIP, entanglement is a 
uniquely quantum phenomenon. Entanglement is a con- 
jectured signature of quantum chaos H, S IM El and 
plays an important role in studies of decoherence |l(J and 
measurement. Understanding entanglement allows us to 
better exploit it as a QIP resource and provides insight 
into the working of quantum mechanics. 

In this paper we study the entanglement production 
of operator classes that approach, but do not properly 
cover, CUE. The purpose of this is two-fold. First, it al- 
lows us to explore the relationship between randomness 
and entanglement and investigate which statistical prop- 
erties of randomness lead to entanglement production. 
Second, some of the operators explored here may be im- 
plemented efficiently. Thus, this study is an exercise of 
how best to produce highly entangled random states. 

The first class of operators we explore are ensembles 
of random matrices which interpolate between integrable 
and CUE 12]. These operators require the same num- 
ber of random parameters as CUE but the parameters 
are drawn from restricted intervals. The second class is 
pseudo-random (PR) operators 0, 0, 0] , possibly effi- 
cient substitutes for random operators in QIP protocols. 
These operators fall far short of the requisite number of 
random parameters when compared to CUE operators. 
The third class is quantum analogs of classically chaotic 
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systems. These operators are the most restricted in terms 
of random parameters but are known to exhibit certain 
statistical properties of random matrices |l5| and can be 
efficiently implemented on a quantum computer. 

As a practical measure of multi-partite entanglement, 
we use the average bipartite entanglement between each 
qubit and the rest of the system 0, , 
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the norm-squared of the wedge product between \v,j) 
and \vj), and pj is the reduced density matrix of qubit 
j. We apply matrices from the above classes to com- 
putational basis states and the average entanglement 
produced as a function of time (number of iterations), 
(Q(t)), is compared to the CUE average entanglement 
(Q)cue = (N — 2)/(N + 1) 18J, where N is the Hilbert 
space dimension. Other entanglement measures, specif- 
ically the concurrence between the two most significant 
qubits and linear entropy between the two N/2 dimen- 
sional subspaces, exhibit behavior similar to Q for the 
operators explored here. 

The statistical properties we examine are the level 
or number variance, the randomness of the eigenvec- 
tors, and the matrix element distribution. The num- 
ber variance measures a two-point eigenvalue correlation 
function which, for many dynamical systems, is known 
to deviate from CUE at long range due to short peri- 
odic orbits [lflj j. Thus, the number variance provides in- 
sight into entanglement generation as a function of time 
since periodic orbits will cause deviations from CUE. In 
the limit of large N, the CUE number variance is [i(J 
Y? CUE {L) = 4^ (In(27rL) + 1 + 7), where 7 ~ .577 is the 
Euler constant. 

A lower bound for the asymptotic bipartite entangle- 
ment production with respect to time, S asy , is the bipar- 
tite entanglement of the systems' eigenvectors, £>eiq mi- 
nus one [2lJ. This result can be extended to Q since it is 
an average of bipartite entanglements, Qasy > 2Qeig ~ 1- 
Let c l k denote the fcth component of the Zth system eigen- 
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vector. The distribution of amplitudes, r\ — c' fc | 2 , for 
CUE eigenvectors in the limit N — > oo and after rescal- 
ing to unit mean is Pcue{v) — e ~ y \ where y — Nrj [2^ . 

When applying an operator to a computational basis 
state the resulting state is a column of the applied op- 
erator. Repeated applications are simply powers of the 
column. Thus, an operator's matrix elements, the result- 
ing state elements, play a central role in the amount of 
entanglement generated. This is seen by writing the av- 
erage Q over all states in terms of |ci| 2 , the state elememt 
amplitudes 

N/2 N N/2 
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CUE matrices can be generated by multiplying eigenvec- 
tors of a Gaussian unitary ensemble (GUE) Hcrmitian 
matrix by random phases and using the resulting vectors 
as matrix columns 23] . Since the eigenvector distribu- 
tion of CUE and GUE are the same 0] and multipli- 
cation by a phase does not change the amplitude of the 
elements, Pcue(%)i the distribution of the rescaled am- 
plitude of CUE matrix elements x, is equal to Pcue{v)- 
The closeness of an operator's matrix element amplitude 
distribution to that of CUE indicates of how much en- 
tanglement the operator can generate. 

The interpolating ensembles are a one-parameter in- 
terpolation between diagonal matrices with uniform, in- 
dependently distributed elements, and CUE 12] . They 
are constructed based on the Hurwitz parameterization 
of CUE matrices. The CUE construction starts with el- 
ementary unitary transformations, E^'^ (<p, ip, x), with 
non-zero elements (23[ |^| 
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which are used to form N — 1 composite rotations 
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and, finally, Ucue = e lot E\E 2 . . .En-i- Angles ip, X, 
and a are drawn uniformly from the intervals 

< tp rs < 2tt, < Xs < 2tt, < a < 2tt, (5) 

and <p rs — sin~ 1 (£ rs 1 ^ 2r+2 - ) ), with £ rs drawn uniformly 

(i i) 

from to 1. The 2x2 block with m, n = i, j and 
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FIG. 1: (Color online) Absolute value of difference be- 
tween the CUE average entanglement, {Q)cue and the av- 
erage Q for interpolating ensemble operators, {Q)s, applied 
t times to each computational basis state, (left) N = 256, 
{Q)cue = .9883 and <5 = .5 (A), .8 (x), .9 (O), -94 (□), 
.96 (v), -98 (o), .99 (left-triangles), .999 (six-pointed stars), 
and .9999 (*). The dotted lines show Qasy determined by the 
eigenvectors of each set of operators. The dashed line shows 
the lower-bound for CUE eigenvectors, (right) N = 256 (□), 
128 (x), 64 (O), and 32 (A). As N decreases the difference 
between {Q)s and (Q)cue as a function of time goes from 
power-law to exponential (for S — .8) and the exponential 
rate increases. This is due to the increased randomness of the 
matrix elements. 

r = is a random SU(2) rotation with respect to the 
Haar measure. Interpolating ensemble construction is the 
same with the angles drawn from constricted intervals 

< tprs < 2tt<S, < Xs < 2tt<5, < a < 2it5, (6) 

with (p rs = sin~ 1 (<5£ r;j 1 ^ 2r+2 ^) and !; rs drawn from to 1. 
The whole is multiplied by a diagonal matrix of random 
phases drawn uniformly from to 2-7T. The parameter S 
ranges from to 1 and provides a smooth transition of 
certain statistical properties between the diagonal circu- 
lar Poisson ensemble and CUE ^] ■ 

For our purposes the interpolating ensembles have the 
same number of random parameters as CUE matrices, 
drawn, however, from restricted intervals. We stress that 
(5 restricts all iV 2 independent variables, such that even 
ensembles of the highest S used here cover only an ex- 
tremely small fraction of CUE space. Figure Q] shows the 
difference between (Q)$, and (Q)cue as a function of 
time and 5. For 5 > .96, (Q)g approaches (Q)cue as an 
exponential whose rate decreases with decreasing S. For 
5 < .5, the approach is a power-law. For constant 5, {Q)s 
approaches its asymptotic value faster for lower N while 
the difference between the asymptotic value and {Q)cue 
increases for lower N. 

The next class of operators we investigate is PR op- 
erators 0, H3I H3] , potentially efficient replacements of 
CUE operators for QIP. To implement a PR operator 
apply to iterations of the n qubit gate: random SU(2) 
rotation, Eqs. © and (JSJ, to each qubit, evolve the sys- 
tem via all nearest neighbor couplings 11]. The coupling 

operator is U nnc = e l ^ /i] £f=i , where a{ is the 

jth qubit z-direction Pauli spin operator. The random 
rotations are different for each qubit and each iteration. 
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FIG. 2: (Color online) Absolute value of difference between 
{Q)cue and (Q(t)) m . (left) n = 8 qubits, m = 2 (+), 4 (x), 
8 (□), 16 (o), 24 (o), 32 (A), and 40 (*). Each average is 
taken over 100 operators applied to all N = 256 computa- 
tional basis states. The dotted lines show Q a s V as determined 
from the eigenvectors. Note that the lower bound for m > 16 
operators is higher than the lower bound of random eigenvec- 
tors, (right) The same for n = 9 (□), 8 (x), 7 (O), and 6 (A). 
The entanglement production compared to CUE changes only 
slightly as a function of n. 



After the roth iteration, a final set of random rotations 
is applied. 

The total number of random parameters used to create 
a PR operator is 3n(m + 1) + 1 where n is the number of 
qubits. This is compared to 2 2 ™ = N 2 random parame- 
ters needed for a CUE matrix. Unless ro is exponential 
in n these operators cannot cover CUE simply because 
there are too few random parameters. 

The absolute value of the difference between (Q)cue 
and (Q) m for n = 8 PR operators, as a function of time 
is shown in Fig. [21 As with the interpolating ensem- 
ble operators, (Q) m approaches (Q)cue as a power-law, 
for low values of ro, less coverage of CUE, and as an 
exponential for greater m. There are a number of inter- 
esting features in this plot. First, for values of ro where 
(Q(t)) m approaches (Q)cue exponentially, the average 
entanglement fluctuates around (Q)cue after the expo- 
nential saturates (Fig. [21plots absolute value). (Q)4o con- 
verges immediately into these fluctuations so increasing 
ro beyond 40 will not increase entanglement generation. 
Also, for operators exhibiting exponential convergence, 
an operator with m = mi at time ti has approximately 
the same (Q) as an operator with m,2 = ami at time 
t-2 = ti/a. For example, (Q(t = 1)) OT =24 is about equal 
to {Q(t — 3)) m= s- This is not the case for operators 
that exhibit non-exponential decay. To create states with 
(Q) ~ (Q)cue one can apply an m = 40 operator once 
or an ro, = 8 operator 5 times. These procedures take the 
same amount of time, but the m — 8 operator requires 
fewer random parameters. Thus, the number of indepen- 
dent variables needed to create entanglement ~ (Q)cue 
in a reasonable time is the number required for the low- 
est m that gives exponential convergence (for n = 8 this 
is approximately m — 8). Unlike the interpolating en- 
semble matrices, the behavior of Q as a function of time 
barely changes with N. 



The above shows that one can create states with CUE 
levels of multi-partite entanglement though only a small 
portion of CUE is covered. An m = 8 PR operator, 
for example, has only 193/256 2 = .3% of the random 
parameters needed for CUE operators but can generate 
CUE levels of entanglement by iterating the operator 5-6 
times. 

The interpolating ensemble operators and PR opera- 
tors lead to similar average entanglement generation be- 
havior as a function of time. For both the entanglement 
approaches (Q)cue as a power-law and, as the operators 
cover more of CUE, an exponential. A priori, there is 
no reason that different restrictions on CUE should give 
rise to similar average entanglement behavior, especially 
since the distributions of Q after one iteration are very 
different for the two types of operators (not shown). 

In light of these results Fig. [3]shows how various statis- 
tical properties relate to entanglement production with 
the aim: to explain why the entanglement generation ap- 
proaches (Q)cue as a power-law or exponential, and why 
(Q)s is dependent on N while (Q)pr is not. 

Based on our numerical investigations, the number 
variance determines the approach of (Q) to (Q)cue- For 
the interpolating ensembles with S < .9 and low ro PR 
operators the number variance diverges almost immedi- 
ately from T, 2 JUE (L). For values of 5 and m where (Q) 
approaches (Q)cue exponentially as a function of time, 
£jj(L) follows CUE faithfully even for large L. In other 
words, operators for which (Q) approaches (Q)cue as a 
power-law do not follow T,q UE (L), while those that ap- 
proach (Q)cue as an exponential follow T, 2 - ;UE (L) up to 
long range correlations, corresponding to short time. 

As expressed in Eq. [21 the matrix element distribution 
(Fig. [2| is key in determining entanglement generation. 
The matrix element distribution approaches the CUE dis- 
tribution at a rate similar to that of the one iteration en- 
tanglement distribution, P(Q) (not shown). This rate of 
convergence is slower than that of other explored statis- 
tics. Second, the matrix element distributions explains 
the behavior of the entanglement generation as a func- 
tion of N. For the interpolating ensemble operators, the 
eigenvector and eigenvalue statistics are practically un- 
changed with N [l2|. The matrix element distribution, 
however, is strongly dependent on N as is the entangle- 
ment generation. For 5 — .9 operators, the matrix ele- 
ment distribution is fully random for N = 8 but deviates 
as N increases. For the PR operators the eigenvector 
and eigenvalue statistics depend only slightly on N, as 
with the matrix element distribution. The entanglement 
generation is also essentially unchanged. 

The final class of operators we explore are quantum 
chaotic operators. While a full analysis is beyond the 
scope of this paper we discuss them in comparison to the 
other operators. Figure[3shows entanglement generation 
as a function of time for the quantum baker's map j2{| 
and an ensemble of chaotic sawtooth [2(|, and Harper 
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FIG. 3: (Color online) (left) Number variance (left) and ma- 
trix element distribution (right), for interpolating ensemble 
(top) and PR operators (bottom). Shown are statistics for 
interpolating ensembles with S = .1 (dashed line), .5 (dot- 
ted line), .9 (chained line) and .98 (light solid line). For the 
PR operators statistics are shown for m — 2 (+), 4 (x), 8 
(□), and 16 (O)- As m and 5 are increased the statistical 
properties approach CUE (solid line). The rate of approach 
for the matrix element distribution is slower than other statis- 
tics, (right) Absolute value of difference between (Q)cue and 
(Q(t)) of the quantum baker's map (o) and sawtooth map for 
non-integer kick strength < k saw < 5 (O) an d Harper's 
map for 1 < "/h < 6 (*). 

maps 27.) . The convergence to (Q)cue is between expo- 
nential and Gaussian but the asymptotic value is lower 
than the other operators. The chaotic maps require on 
order n 2 [25| or n 3 gates. The PR operators re- 
quire n — 1 coupling terms and n rotations per iteration, 
approximately 2mn gates per operator. These are com- 
parable as long as m is less than quadratic in n. How- 
ever, unlike the quantum chaotic operators, the coupling 
terms at every iteration of the PR algorithm can be ap- 
plied simultaneously. The PR operators thus require only 
2m + 1 gates which appears to be less then that required 
of chaotic operators. 

In conclusion, we have studied the entanglement gen- 
eration of interpolating ensemble, PR, and quantum 
chaotic operators as a function of time. These opera- 
tors restrict the full space of CUE in different ways and 
the effect of these restrictions can be seen in the entan- 
glement generation. The statistical properties which in- 
fluence the entanglement generation include the number 
variance, which effects the entanglement generation as 
a function of time, and the matrix element distribution, 
which determines other aspects of the entanglement. Fi- 
nally, we note that the PR operators may be efficiently 
implemented on a quantum computer and provide a way 
to create highly entangled, random states. 
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